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Abstract 

This paper provides a summary of the current equations and rescaling factors for converting 
calibrated Digital Numbers (DNs) to absolute units of at-sensor spectral radiance, Top-Of- 
Atmosphere (TOA) reflectance, and at-sensor brightness temperature. It tabulates the necessary 
constants for the Multispectral Scanner (MSS), Thematic Mapper (TM), Enhanced Thematic 
Mapper Pius (ETM+), and Advanced Land imager (ALI) sensors. These conversions provide a 
basis for standardized comparison of data in a single scene or between images acquired on 
different dates or by different sensors. This paper forms a needed guide for Landsat data users 
who now have access to the entire Landsat archive at no cost. 
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1. Introduction 

The Landsat series of satellites provides the longest continuous record of satellite-based 
observations. As such, Landsat is an invaluable resource for monitoring global change and is a 
primary source of medium spatial resolution Earth observations used in decision-making (Fuller et 
al., 1994; Townshend et al., 1995; Goward et al., 1997; Vogelmann et al., 2001; Woodcock et al., 
2001; Cohen et al., 2004; Goward et al., 2006; Masek et al., 2008; Wulder et al., 2008). To meet 
observation requirements at a scale revealing both natural and human-induced landscape 
changes, Landsat provides the only inventory of the global land surface over time on a seasonal 
basis (Special issues on Landsat, 1984; 1985; 1997; 2001; 2003; 2004; 2006). The Landsat 
Program began in early 1972 with the launch of the first satellite in the series. As technological 
capabilities increased, so did the amount and quality of image data captured by the various 
sensors onboard the satellites. Table 1 presents general information about each Landsat satellite. 

Landsat satellites can be classified into three groups, based on sensor and platform 
characteristics. The first group consists of Landsat 1 (LI), Landsat 2 (L2), and Landsat 3 (L3), 
with the Multispectral Scanner (MSS) sensor and the Return Beam Vidicon (RBV) camera as 
payloads on a “NIMBUS-like” platform. The spatial resolution of the MSS sensor was 
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approximately 79 m (but often processed to pixel size of 60 m), with four bands ranging from the 
visible blue to the Near-Infrared (NIR) wavelengths. The MSS sensor on L3 included a fifth band 
in the thermal infrared wavelength, with a spectral range from 10.4 to 12.6 pm. The L1-L3 MSS 
sensors used a band-naming convention of MSS-4, MSS-5, MSS-6, and MSS-7 for the blue, 
green, red, and NIR bands, respectively (Markham & Barker, 1983). This designation is obsolete, 
and to be consistent with the TM and ETM+ sensors, the MSS bands are referred to here as 
Bands 1-4, respectively. 

The second group includes Landsat 4 (L4) and Landsat 5 (L5), which carry the Thematic 
Mapper (TM) sensor, as well as the MSS, on the Multimission Modular Spacecraft. This second 
generation of Landsat satellites marked a significant advance in remote sensing through the 
addition of a more sophisticated sensor, improved acquisition and transmission of data, and more 
rapid data processing at a highly automated processing facility. The MSS sensor was included to 
provide continuity with the earlier Landsat missions, but TM data quickly became the primary 
source of information used from these satellites because the data offered enhanced spatial, 
spectral, radiometric, and geometric performance over data from the MSS sensor. The TM sensor 
has a spatial resolution of 30 m for the six reflective bands and 120 m for the thermal band. 
Because there are no onboard recorders on these sensors, acquisitions are limited to real-time 
downlink only. 

The third group consists of Landsat 6 (L6) and Landsat 7 (L7), which include the 
Enhanced Thematic Mapper (ETM) and the Enhanced Thematic Mapper Plus (ETM+) sensors, 
respectively. No MSS sensors were included on either satellite. Landsat 6 failed on launch. The 
L7 ETM+ sensor has a spatial resolution of 30 m for the six reflective bands, 60 m for the thermal 
band, and includes a panchromatic (pan) band with a 15 m resolution. L7 has a 378 gigabit (Gb) 
Solid State Recorder (SSR) that can hold 42 minutes (approximately 100 scenes) of sensor data 
and 29 hours of housekeeping telemetry concurrently (L7 Science Data User's Handbook 1 ). 

The Advanced Land Imager (ALI) onboard the Earth Observer-1 (EO-1) satellite is a 
technology demonstration that serves as a prototype for the Landsat Data Continuity Mission 
(LDCM). The ALI observes the Earth in 10 spectral bands; nine spectral bands have a spatial 
resolution of 30 m, and a pan band has a spatial resolution of 10 m. 

The Landsat data archive at the U.S. Geological Survey (USGS) Earth Resources 
Observation and Science (EROS) Center holds an unequaled 36-year record of the Earth’s 


1 http 1 v'/landsathancibo ok.qsfc.nasa qoy/handbook.htmj . Landsat Project Science Office, Goddard Space Flight Center. 
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surface and is available at no cost to users via the Internet (Woodcock et al., 2008). Users can 
access and search the Landsat data archive via the EarthExplorer (EE) 2 or Global Visualization 
Viewer (GloVis) 3 web sites. Note that the Landsat scenes collected by locations within the 
International Ground Station (IGS) network may be available only from the particular station that 
collected the scene. 

2. Purpose 

Equations and parameters to convert calibrated Digital Numbers (DNs) to physical units, 
such as at-sensor radiance or Top-Of-Atmosphere (TOA) reflectance, have been presented in a 
“sensor-specific” manner elsewhere, e.g., MSS (Markham & Barker, 1986, 1987; Helder, 1993), 
TM (Chander & Markham, 2003; Chander et al., 2007), ETM+ (Handbook 1 ), and ALI (Markham et 
al., 2004a). This paper, however, tabulates the necessary constants for all of the Landsat sensors 
in one place defined in a consistent manner and provides a brief overview of the radiometric 
calibration procedure summarizing the current accuracy of the at-sensor spectral radiances 
obtained after performing these radiometric conversions on standard data products generated by 
U S. ground processing systems. 

3. Radiometric calibration procedure 

The ability to detect and quantify changes in the Earth's environment depends on 
sensors that can provide calibrated (known accuracy and precision) and consistent 
measurements of the Earth's surface features through time. The correct interpretation of scientific 
information from a global, long-term series of remote-sensing products requires the ability to 
discriminate between product artifacts and changes in the Earth processes being monitored (Roy 
et al., 2002). Radiometric characterization and calibration is a prerequisite for creating high- 
quality science data, and consequently, higher-level downstream products. 

3.1. MSS sensors 

Each MSS sensor incorporates an Internal Calibrator (1C) system, consisting of a pair of 
lamp assemblies (for redundancy) and a rotating shutter wheel. The shutter wheel includes a 
mirror and a neutral density filter that varies in transmittance with rotation angle. The calibration 
system output appears as a light pulse at the focal plane that rises rapidly and then decays 
slowly. This pulse is referred to as the calibration wedge (Markham & Barker, 1987). The 
radiometric calibration of the MSS sensors is performed in two stages. First, raw data from Bands 
1-3 are "decompressed" or linearized and rescaled to 7 bits using fixed look-up tables. The look- 


2 http - '/ear t hexpiorer.usqs.qov 

3 h!t o:/'Qloyis usgs.go v 
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up tables are derived from prelaunch measurements of the compression amplifiers. Second, the 
postlaunch gain and offset for each detector of all four bands are individually calculated by a 
linear regression of the detector responses to the samples of the in-orbit calibration wedge with 
the prelaunch radiances for these samples. A reasonable estimate of the overall calibration 
uncertainty of each MSS sensor at-sensor spectral radiances is ±10%, which was the specified 
accuracy for the sensor (Markham & Barker, 1987). In most cases, the ground processing system 
must apply an additional step to uncalibrate the MSS data because a number of MSS scenes 
were archived as radiometrically corrected products. The previously calibrated archived MSS 
data must be transformed back into raw DNs using the coefficients stored in the data before 
applying the radiometric calibration procedure. Studies are underway to evaluate the MSS 
calibration consistency and provide post-calibration adjustments of the MSS sensors so they are 
consistent overtime and consistent between sensors (Helder, 2008a). 

3.2. TM sensors 

The TM sensor includes an onboard calibration system called the 1C. The 1C consists of a 
black shutter flag, three lamps, a cavity blackbody, and the optical components necessary to get 
the lamp and blackbody radiance to the focal plane. The lamps are used to calibrate the reflective 
bands, and the blackbody is used to calibrate the thermal band. Historically, the TM radiometric 
calibration procedure used the detector’s response to the 1C to determine radiometric gains and 
offsets on a scene-by-scene basis. Before launch, the effective radiance of each lamp state for 
each reflective band’s detector was determined such that each detector’s response to the internal 
lamp was compared to its response to an external calibrated source. The reflective band 
calibration algorithm for in-flight data used a regression of the detector responses against the 
prelaunch radiances of the eight lamp states. The slope of the regression represented the gain, 
while the intercept represented the bias. This algorithm assumed that irradiance of the calibration 
lamps remained constant over time since launch. Any change in response was treated as a 
change in sensor response, and thus was compensated for during processing. On-orbit data from 
individual lamps indicated that the lamps were not particularly stable. Because there was no way 
to validate the lamp radiances once in orbit, the prelaunch measured radiances were the only 
metrics available for the regression procedure, Recent studies 4 (Thome et al., 1997a, 1997b; 
Helder et al., 1998; Markham et al., 1998; Teillet et al., 2001, 2004; Chander et al., 2004) indicate 
that the regression calibration did not actually represent detector gains for most of the mission. 
However, the regression procedure was used until 2003 to generate L5 TM data products and is 
still used to generate L4 TM products. The calibration uncertainties of the L4 TM at-sensor 


4 

Radiometric performance studies of the TM sensors have also led to a detailed understanding of several image artifacts 
due to particular sensor characteristics (Helder & Ruggles, 2004a). These artifact corrections (such as Scan-Correlated 
Shift [SCS], Memory Effect [ME], and Coherent Noise [CN]), along with detector-to-detector normalization (Helder et. al., 
2004b), are necessary to maintain the internal consistency of the calibration within a scene. 
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spectral radiances are ±10%, which was the specified accuracy for the sensor (GSFC 
specification, 1981). 

The L5 TM reflective band calibration procedure was updated in 2003 (Chander & 
Markham, 2003) to remove the dependence on the changing 1C lamps. The new calibration gains 
implemented on May 5, 2003, for the reflective bands (1-5, 7) were based on lifetime radiometric 
calibration curves derived from the detectors' responses to the 1C, cross-calibration with ETM+, 
and vicarious measurements (Chander et al., 2004a). The gains were further revised on April 2, 
2007, based on the detectors’ responses to pseudo-invariant desert sites and cross-calibration 
with ETM+ (Chander et al., 2007). Although this calibration update applies to all archived and 
future L5 TM data, the principal improvements in the calibration are for data acquired during the 
first eight years of the mission (1984-1991), where changes in the sensor gain values are as 
much as 15%. The radiometric scaling coefficients for Bands 1 and 2 for approximately the first 
eight years of the mission have also been changed. Along with the revised reflective band 
radiometric calibration on April 2, 2007, an sensor offset correction of 0.092 W/(m 2 sr pm), or 
about 0.68 K (at 300 K), was added to all L5 TM thermal band (Band 6) data acquired since April 
1999 (Barsi et al., 2007). The L5 TM radiometric calibration uncertainty of the at-sensor spectral 
radiances is around 5% and is somewhat worse for early years, when the sensor was changing 
more rapidly, and better for later years (Helder et al., 2008b). The L4 TM reflective bands and the 
thermal band on both the TM sensors continue to be calibrated using the 1C. Further updates to 
improve the thermal band calibration are being investigated, as is the calibration of the L4 TM. 

3.3. ETM+ sensor 

The ETM+ sensor has three onboard calibration devices for the reflective bands: a Full 
Aperture Solar Calibrator (FASC), which is a white painted diffuser panel; a Partial Aperture Solar 
Calibrator (PASC), which is a set of optics that allows the ETM+ to image the Sun through small 
holes; and an 1C, which consists of two lamps, a blackbody, a shutter, and optics to transfer the 
energy from the calibration sources to the focal plane. The ETM+ sensor has also been calibrated 
vicariously using Earth targets such as Railroad Valley (Thome, 2001; Thome et al., 2004) and 
cross-calibrated with multiple sensors (Teillet et al,, 2001, 2006, 2007; Thome et al., 2003; 
Chander et al., 2004b, 2007b, 2008). The gain trends from the ETM+ sensor are regularly 
monitored on-orbit using the onboard calibrators and vicarious calibration. The calibration 
uncertainties of ETM+ at-sensor spectral radiances are ±5%. ETM+ is the most stable of the 
Landsat sensors, changing by no more than 0.5% per year in its radiometric calibration (Markham 
et al., 2004b). The ETM+ radiometric calibration procedure uses prelaunch gain coefficients 
populated in the Calibration Parameter File (CPF). These CPFs, issued quarterly, have both an 
“effective" and “version” date. The effective date of the CPF must match the acquisition date of 
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the scene. A CPF version is active until a new CPF for that date period supersedes it. Data can 
be processed with any version of a CPF; the later versions have more refined parameters, as 
they reflect more data-rich post-acquisition analysis. 

The ETM+ images are acquired in either a low- or high-gain state. The goal of using two 
gain settings is to maximize the sensors’ 8-bit radiometric resolution without saturating the 
detectors. For all bands, the low-gain dynamic range is approximately 1 .5 times the high-gain 
dynamic range. Therefore, low-gain mode is used to image surfaces with high brightness (higher 
dynamic range but low sensitivity), and high-gain mode is used to image surfaces with low 
brightness (lower dynamic range but high sensitivity). 

All of the ETM+ acquisitions after May 31, 2003, have an anomaly caused by the failure 
of the Scan Line Corrector (SLC), which compensated for the forward motion of the spacecraft so 
that all the scans were aligned parallel with each other. The images with data loss are referred to 
as SLC-off images, whereas images collected prior to the SLC failure are referred to as SLC-on 
images (i.e. , no data gaps exist). The malfunction of the SLC mirror assembly resulted in the loss 
of approximately 22% of the normal scene area (Storey et al„ 2005). The missing data affects 
most of the image, with scan gaps varying in width from one pixel or less near the center of the 
image to 14 pixels along the east and west edges of the image, creating a repeating wedge- 
shaped pattern along the edges. The middle of the scene, approximately 22 km wide on a Level 1 
product, contains very little duplication or data loss. Note that the SLC failure has no impact on 
the radiometric performance with the valid pixels. 

3.4. ALI sensor 

The ALI has two onboard radiometric calibration devices: a lamp-based system and a 
solar-diffuser with variable irradiance controlled by an aperture door. In addition to its onboard 
calibrators, ALI has the ability to collect lunar and stellar observations for calibration purposes. 
The ALI radiometric calibration procedure uses a fixed set of detector-by-detector gains 
established shortly after launch and biases measured shortly after each scene acquisition by 
closing the ALI’s shutter. The calibration uncertainties of the ALI at-sensor spectral radiances are 
±5% (Mendenhall & Lencioni, 2002). The ALI sensor is well-behaved and stable, with changes in 
the response being less than 2% per year even early in the mission, and averaging, at most, 
slightly more than 1% per year over the full mission (Markham et al., 2006). 

4. Conversion to at-sensor spectral radiance (Q car to- Lj.) 

Calculation of at-sensor spectral radiance is the fundamental step in converting image 
data from multiple sensors and platforms into a physically meaningful common radiometric scale. 
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Radiometric calibration of the MSS, TM, ETM+, and ALI sensors involves rescaling the raw digital 
numbers (Q) transmitted from the satellite to calibrated digital numbers (Q ca i) 5 , which have the 
same radiometric scaling for all scenes processed on the ground for a specific period. 

During radiometric calibration, pixel values (Q) from raw, unprocessed image data are 
converted to units of absolute spectral radiance using 32-bit floating-point calculations. The 
absolute radiance values are then scaled to 7-bit (MSS, Q ca imax = 127), 8-bit (TM and ETM+, 
Qcaimax = 255), and 16-bit (ALI, Q ca i max = 32767) numbers representing Q cat before output to 
distribution media. Conversion from Q ca i in Level 1 products back to at-sensor spectral radiance 
(L x ) requires knowledge of the lower and upper limit of the original rescaling factors. The following 
equation is used to perform the Q ca i-to- L, conversion for Level 1 products: 


= 


or 
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Where 

Lx = Spectral radiance at the sensor’s aperture [W/(m 2 sr pm)] 

Q ca i = Quantized calibrated pixel value [DN] 

Qcaimin = Minimum quantized calibrated pixel value corresponding to LMIN,[DN] 
Qcaimax = Maximum quantized calibrated pixel value corresponding to LMAX, [DN] 
LMIN, = Spectral at-sensor radiance that is scaled to Q ca imin [W/(m 2 sr pm)] 

LMAX, = Spectral at-sensor radiance that is scaled to Q ca imax [W/(m 2 sr pm)] 

Q rescale = Band-specific rescaling gain factor [(W/(m 2 sr pm))/ DN] 

Brescaie = Band-specific rescaling bias factor [W/(m 2 sr pm)] 


( 1 ) 


Historically, the MSS and TM calibration information is presented in spectral radiance 
units of mW/(cm 2 sr pm). To maintain consistency with ETM+ spectral radiance, units of W/(m z sr 
pm) are now used for MSS and TM calibration information. The conversion factor is 1:10 when 

5 These are the DNs that users receive with Level 1 Landsat products. 
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converting from mW / (cm 2 sr pm) units to W/ (m 2 sr pm). Tables 2, 3, 4, and 5 summarize the 
spectral range, post-calibration dynamic ranges 6 (LMIN A and LMAX ; . scaling parameters and the 
corresponding rescaling gain [G reS caie]] and rescaling bias [B res cate] values), and mean 
exoatmospheric solar irradiance (ESUN A ) for the MSS, TM, ETM+, and ALI sensors, respectively. 

Tables 2-5 give the prelaunch “measured” (as-built performance) spectral ranges. These 
numbers are slightly different from the original filter specification. The center wavelengths are the 
average of the two spectral range numbers. Figures 1 and 2 show the Relative Spectral 
Response (RSR) profiles of the Landsat MSS (Markham & Barker, 1983), TM (Markham & 
Barker, 1986), ETM+ (Handbook 1 ), and ALI (Mendenhall & Parker, 1999) sensors measured 
during prelaunch characterization. The ETM+ spectral bands were designed to mimic the 
standard TM spectral bands 1-7. The ALI bands were designed to mimic the six standard ETM+ 
solar reflective spectral bands 1-5, and 7; three new bands, Ip, 4p, and 5p, were added to more 
effectively address atmospheric interference effects and specific applications. The ALI band 
numbering corresponds with the ETM+ spectral bands. Bands not present on the ETM+ sensor 
are given the “p,” or prime, designation. MSS spectral bands are significantly different from TM 
and ETM+ spectral bands. 

The post-calibration dynamic ranges are band-specific rescaling factors typically provided 
in the Level 1 product header file. Over the life of the Landsat sensors, occasional changes have 
occurred in the post-calibration dynamic range. Future changes are anticipated, especially in the 
MSS and TM data, because of the possible adjustment of the calibration constants based on 
comparisons to absolute radiometric measurements made on the ground. In some cases, the 
header file may have different rescaling factors than provided in the table included here. In these 
cases, the user should use the header file information that comes with the product. 

Two processing systems will continue to generate Landsat data products: the Level 1 
Product Generation System (LPGS) and the National Land Archive Production System (NLAPS). 
Starting December 8, 2008, all L7 ETM+ and L5 TM (except Thematic Mapper-Archive [TM-A] 7 
products) standard Level 1 products are processed through the LPGS, and all L4 TM and MSS 


6 The post-calibration dynamic ranges summarized in Tables 2-5 are only applicable to Landsat data processed and 
distributed by the USGS EROS Center. The IGSs may process the data differently, and these rescaling factors may not 
be applicable. “Special collections," such as the Multi-Resolution Land Characteristics Consortium (MRLC) or Global Land 
Survey (GLS), may have a different processing history, so the user needs to verify the respective product header 
information. 

' A small number of TM scenes were archived as radiometricalty corrected products known as TM-A data. The TM-A data 
are archived on a scene-by-scene basis (instead of intervals). The L4 and L5 TM-A scenes will continue to be processed 
using NLAPS (with Q ca i n ,„ =0), which attempts to uncalibrate the previously applied calibration and generates the product 
using updated calibration procedures. Note that approximately 80 L4 TM and approximately 13,300 L5 TM scenes are 
archived as TM-A data, with acquisition dates ranging between Sept.1982 and Aug. 1990. 
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standard Level 1 products are processed through the NLAPS. The Landsat Program is working 
toward transitioning the processing of all Landsat data to LPGS (Kline, personal communication). 
In mid-2009, the processing of L4 TM data will transition from NLAPS to LPGS. The scenes 
processed using LPGS include a header file (.MTL), which lists the LMIN A and LMAX>. values but 
not the rescaling gain and bias numbers. The scenes processed using NLAPS include a 
processing history work order report (.WO), which lists the rescaling gain and bias numbers but 
not the LMIN X and LMAX x . 

The sensitivity of the detector changes over time, causing a change in the detector gain 
applied during radiometric calibration. However, the numbers presented in Tables 2-5 are the 
rescaling factors, which are the post-calibration dynamic ranges. The LMIN* and LMAX a are a 
representation of how the output Landsat Level 1 data products are scaled in at-sensor radiance 
units. Generally, there is no need to change the LMIN* or LMAX* unless something changes 
drastically on the sensor. Thus, there is no time dependence for any of the rescaling factors in 
Tables 2-5. 

5. Conversion to TOA reflectance (L x -to- p P ) 

A reduction in scene-to-scene variability can be achieved by converting the at-sensor 
spectral radiance to exoatmospheric TOA reflectance, also known as in-band planetary albedo. 
When comparing images from different sensors, there are three advantages to using TOA 
reflectance instead of at-sensor spectral radiance. First, it removes the cosine effect of different 
solar zenith angles due to the time difference between data acquisitions. Second, TOA 
reflectance compensates for different values of the exoatmospheric solar irradiance arising from 
spectral band differences. Third, the TOA reflectance corrects for the variation in the Earth-Sun 
distance between different data acquisition dates. These variations can be significant 
geographically and temporally. The TOA reflectance of the Earth is computed according to the 
equation: 

k • L, d 2 

p, = (2) 

ESUN, • cos Q s 

where 

p x = Planetary TOA reflectance [unitless] 

tt = Mathematical constant approximately equal to 3.14159 [unitless] 

L x =Spectral radiance at the sensor’s aperture [W/(m 2 sr pm)] 
d = Earth-Sun distance [astronomical units] 

ESUN* = Mean exoatmospheric solar irradiance [W/(m 2 pm)] 
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0 S = Solar zenith angle [degrees 8 ] 

Note that the cosine of the solar zenith angle is equal to the sine of the solar elevation 
angle. The solar elevation angle at the Landsat scene center is typically stored in the Level 1 
product header file (.MTL or .WO) or retrieved from the USGS EarthExplorer or GloVis online 
interfaces under the respective scene metadata (these web sites also contain the acquisition time 
in hours, minutes, and seconds). The reflectance calculation requires the Earth-Sun distance (d). 
Table 6 presents d in astronomical units throughout a year generated using the Jet Propulsion 
Laboratory (JPL) Ephemeris 9 10 (DE405) data. The d numbers are also tabulated in the Nautical 
Almanac. 

The last column of Tables 2-5 summarizes solar exoatmospheric spectral irradiances 
(ESUNx) for the MSS, TM, ETM+, and ALI sensors using the Thuillier solar spectrum (Thuillier et 
al., 2003). The Committee on Earth Observation Satellites (CEOS) Working Group on Calibration 
and Validation (WGCV) recommends 0 using this spectrum for applications in optical-based Earth 
Observation that use an exoatmospheric solar irradiance spectrum. The Thuillier spectrum is 
believed to be the most accurate and an improvement over the other solar spectrum. Note that 
the CHKUR solar spectrum in MODTRAN 4.0 (Air Force Laboratory, 1998) was used previously 
for ETM+ (Handbook 1 ) and TM (Chander & Markham, 2003), whereas the Neckel and Lab 
(Neckel & Lab, 1984) and Iqbal (Iqbal, 1983) solar spectrums were used for MSS and TM solar 
irradiance values (Markham & Barker, 1986). The primary differences occur in Bands 5 and 7. 
For comparisons to other sensors, users need to verify that the same solar spectrum is used for 
all sensors. 

6. Conversion to at-sensor brightness temperature (L A -to- T) 

The thermal band data (Band 6 on TM and ETM+) can be converted from at-sensor 
spectral radiance to effective at-sensor brightness temperature. The at-sensor brightness 
temperature assumes that the Earth's surface is a black body (i.e., spectral emissivity is 1), and 
includes atmospheric effects (absorption and emissions along path). The at-sensor temperature 
uses the prelaunch calibration constants given in Table 7. The conversion formula from the at- 
sensor’s spectral radiance to at-sensor brightness temperature is: 


8 Note that Excel, Matlab, C, and many other software applications use radians, not degrees, to perform calculations. The 
conversion from degrees to radians is a multiplication factor of pi/180. 

9 http:/.ssd.ipl.nasa.qov/?horizons 

10 CEOS-recommended solar irradiance spectrum, http 'wgcv ceos orq 
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In 

where: 

T = Effective at-sensor brightness temperature [K] 

K2 = Calibration constant 2 [KJ 

K1 = Calibration constant 1 [W/(m 2 sr pm)] 

l_ A = Spectral radiance at the sensor’s aperture [W/(m 2 sr pm)] 

In = Natural logarithm 

The ETM+ Level 1 product has two thermal bands, one acquired using a low gain setting 
(often referred to as Band 6L; useful temperature range of 130-350 K) and the other using a high 
gain setting (often referred to as Band 6H; useful temperature range of 240-320 K). The noise 
equivalent change in temperature (NEAT) at 280 K for ETM+ high gain is 0.22 and for low gain is 
0.28. The TM Level 1 product has only one thermal band (there is no gain setting on the TM 
sensor), and the thermal band images have a useful temperature range of 200-340 K. The NEAT 
at 280 K for L5 TM is 0.17-0.30 (Barsi et al„ 2003). 

7. Conclusion 

This paper provides equations and rescaling factors for converting Landsat calibrated 
DNs to absolute units of at-sensor spectral radiance, TOA reflectance, and at-sensor brightness 
temperature. It tabulates the necessary constants for the MSS, TM, ETM+, and ALI sensors in a 
coherent manner using the same units and definitions. This paper forms a needed guide for 
Landsat data users who now have access to the entire Landsat archive at no cost. Studies are 
ongoing to evaluate the MSS calibration consistency and provide post-calibration adjustments of 
the MSS sensors so they are consistent over time and consistent between sensors. Further 
updates to improve the TM and ETM+ thermal band calibration are being investigated, as is the 
calibration of the L4 TM. 
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Table 1 

Landsat satellites launch dates 



Sensors 

Launch Date 

Decommission 

[Min 


liianinrt 


Md ?)T”I I'fM 

I 




km 

I decrees ; 



miEwnin 

1 Hit! EMI I 



January 7, 1978 

920 

BH 


18 

9:30 


5 

! 

I 

January 22, 1975 


920 


■nfejfcj!| 

IB 

9:30 

Landsat 3 

MSS and RBV 

March 5, 1978 

■ammt-vai 

920 



18 


Landsat 4 

MSS and TM 

July 16, 1982 


705 

K:KTK‘ 

fefciw 

IB 


Landsat 5 

MSS and TM 

March 1 , 1984 


705 



16 


Landsat 6 

| ETM 

October 5, 1993 

Did not achieve orbit 

1 : n 

Landsat 7 

ETM+ 

April 15, 1999 

Operational 

705 

98.20 


16 

filiT 

EO-1 

ALI 

November 21 , 2000 

Operational 

705 

98.20 

B-Eii 

16 

10 01 1 


Table 2 

MSS spectral range, post-calibration dynamic ranges, and mean exoatmospheric solar irradiance 
(ESUN a ) 


MSS Sensors (G caJmil , - 0 and Q M i max “ 127; 


Band 

Spectral 

Range 

Center 

Wavelength 

LMIN* 

LMAX* 

^rescale 

Brescale 

ESUN* 

Units 

pn» 

| W (m 2 si jim) 

(W in 2 si inn) DN 

W tin 2 si inn) 

W (in 2 . inn) 


LI MSS (NLAPS) 

1 

0.499-0.597 

0.548 

0 

248 

1 .952760 

0 


2 

0.603-0.701 

0.652 

0 

200 

1 .574800 

0 


3 

0 694 -0 800 

0.747 

0 

176 

1.385830 

0 

1276 

4 

0.810-0.989 

0.900 

0 

153 

1 .204720 

0 

880.1 


L2 MSS (NLAPS) 

1 

0.497-0.598 

0.548 

8 

263 

2.007870 

8 

1829 

2 

0.607-0.710 

0.659 

6 

176 

1 .333580 

6 

1539 

3 

0.697 - 0 802 

0.750 

6 

152 

1.149610 

6 

1268 

4 

0 807-0.990 

0 899 

3.66667 

130.333 

0.997373 

3.66667 

886.6 


L3 MSS (NLAPS) 

1 

0.497 -0.593 

WEMSsMM 

4 

259 

2.007870 

4 

1839 

2 

0.606 -0.705 


3 

179 

1.385830 

3 

1555 

3 

0.693 - 0 793 

0.743 

3 

149 

1.149610 

3 

1291 

4 

0.812-0.979 

0.896 

1 

128 

1.000000 

1 

887.9 


L4 MSS (NLAPS) 

1 

0 495 - 0.605 

0.550 

4 

238 


4 


2 

0.603-0.696 

0.650 

4 

164 

1.259840 

4 

■sa 

3 

0.701 -0.813 


5 
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1 .078740 

5 

1260 I 

4 

0 808- 1.023 


4 
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0.881890 

4 

EMBI 


L5 MSS (NLAPS) 

1 

0.497-0 607 

0.552 

3 

268 

2.086610 

3 

1824 

2 

0.603 - 0.697 

0.650 

3 

179 

1 .385830 

3 

1570 

3 

0.704-0.814 

0.759 


warn 

1.125980 

5 

1249 

4 

0.809- 1.036 

0.923 


■EH 

0.944882 


H£X» 


Note 1: In some cases, the header file may have different rescaling factors than provided here. In these cases, the user should 
use the header file information that comes with the product. Tables 1 (Markham 8 Barker, 1986, 1987) provide a summary of the 
band-specific LMIN , and LMAX , rescaling factors that have been used at different times and by different systems for the ground 
processing of MSS data. 
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Table 3 

TM spectral range, post-calibration dynamic ranges, and mean exoatmospheric solar irradiance 
(ESUN a ) 


TM Sensors (Q calmin = 1 and Q calmax = 255) 


Band 

Spectial 

Ramie 

Center 

Wavelength 

LMIN* 


Qrescale 

^rescale 

ESUN h 

Units 

1 l ,m 


(Win 3 si midi DN 

W (m 2 si pm( 

W (m 2 41 inf 

L4 TM (NLAPS) 1 

1 

0.452-0.518 

0.485 


152.10 

0.602431 

-1.52 

1983 

2 

0.529 -0.609 

0.569 

-2.84 

296.81 

1.175098 

-2.84 

1795 

3 

0.624-0.693 

0.659 

-1.17 


0.805765 

-1.17 

1539 

4 

0.776 - 0.905 

0.841 

-1.51 

206.20 

0.814549 

-1.51 

1028 

5 

1 .568 - 1 .784 

1 .676 

-0.37 

27.19 

0.108078 

-0.37 

219.8 

6 

10.42-11.66 

11.040 

1 .2378 

15.3032 

0.055158 

1.2378 

N/A 

7 

2.097-2.347 

2.222 

-0.15 

14.38 

0.056980 

-0.15 

83.49 

L4 TM 

(LPGS) 

1 

0.452-0.518 

0.485 

-1.52 

163 

0.647717 

-2.17 

1983 

-1.52 

171 

0.679213 

-2.20 

2 

0.529-0.609 

0.569 

-2.84 

336 

1.334016 

-4.17 

1795 

3 

0.624 - 0 693 

0.659 

-1.17 

254 

1 .004606 

-2.17 

1539 

4 

0.776-0 905 

0 841 

-1.51 

221 

0.876024 

-2.39 

1028 

5 

1.568-1.784 

1.676 

-0.37 

31.4 

0.125079 

-0.50 

219.8 

6 

10.42-11.66 

11.040 

1 .2378 

15.3032 

0.055376 

1 .2378 

N/A 

7 

2.097-2.347 

2.222 

-0.15 

16.6 

0.065945 

-0 22 

83.49 

L5TM 

(LPGS) i 

1 

0.452-0.518 

0.485 

-1.52 

169 

0.671339 

-219 | 

1983 

-1.52 

193 

0.765827 

-2.29 

2 

0.528 -0.609 

0.569 

-284 

333 

1.322205 


1796 

-2.84 

365 


-4.29 

3 

0.626-0.693 

0.660 

-1.17 

264 

1.043976 

-2.21 

1536 

4 

0.776-0.904 

0.840 

-1.51 

221 


-2.39 

1031 

5 

1.567-1 784 

1.676 

-0.37 

30.2 

0.120354 

-0.49 

220.0 

6 

10.45 - 12.42 

11.435 

1 .2378 

15.3032 

0.055376 

1.18 

N/A 

7 

2.097-2.349 

2.223 

-0.15 

16.5 

0.065551 

-0 22 

83.44 


Note 1: The Q 5 ar»in = 0 for data processed using NLAPS . The Q ca „, n = 1 for data processed using LPGS. 


Note 2: The LMIN , is typically set to a small negative number , so a "zero radiance" target will be scaled to a small positive DN value , 
even In the presence of sensor noise (typically 1 DN or less [1 sigma]). This value Is usually not changed throughout the mission. 
Note 3: In mid-2009, the processing of L4 TM data will transition from NLAPS to LPGS. NLAPS used IC-based calibration. The L4 TM 
data processed by LPGS will be radiometrically calibrated using a new lifetime gain model procedure and revised calibration 
parameters . Use the header file Information that comes with the product and the above rescaling factors will not be applicable . The 
numbers highlighted in grey are the revised (LMAX, = 163.0) post-calibration dynamic ranges for L4 TM Band 1 data acquired 
between July 16,1 982 (launch), and August 23, 1 986. 

Note 4: The radiometric scaling coefficients for L5 TM Bands 1 and 2 for approximately the first eight years (1 984-1 991 ) of the 
mission were changed to optimize the dynamic range and better preserve the sensitivity of the early mission data. The numbers 
highlighted in grey are the revised (LMAX, = 169.0, 333.0) post-calibration dynamic ranges for L5 TM Band 1 8 2 data acquired 
between March 1 , 1 984 (launch), and December 31 , 1 991 (Chander et al., 2007a). 
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Table 4 

ETIVI+ spectral range, post-calibration dynamic ranges, and mean exoatmospheric solar 
irradiance (ESUN A ) 


L7 ETM+ Sensor (Qcaimin — 1 snd Q ca imax = 255) 

Low Gain (LPGS) 1 

Band 

ajaQjgn 

HEU 

Center 

Wavelength 

LMIN„ 

LMAXji 

^rescale 

Brescale 

ESUNj 

Units 

! 

W (m 2 

sr |im) 

(W in 2 sr inn) DN 


MfiMm 

1 

0 . 452 - 0.514 

0.483 

- 6.2 

■55 M 

1.180709 

- 7.38 


2 

0 . 519 - 0.601 

0.560 

kOI 

■c'tifcl 

1 .209843 

- 7.61 


3 

0.631 - 0.692 

0662 

- 5.0 

234.4 

0.942520 

- 5.94 

1533 

4 

0 772 - 0 898 

0.835 

- 5.1 

241.1 

0.969291 

- 6.07 

1039 

5 

1 . 547 - 1.748 

1.648 

- 1.0 

47.57 

0.191220 

- 1.19 

230.8 

G 

10.31 -12 36 

1 1 .335 

00 

17.04 

0.067087 

- 0.07 

N/A 


2.065 - 2.346 

2.206 

- 0.35 

16.54 

0.066496 


84.90 


0 . 515 - 0.896 

0.706 

warn 

243.1 

0.975591 

- 5.68 

1362 



High Gain (LPGS) 


Band 


Center 

Wavelength 

LMINj 

LMAXft 


Brescal© 

ESUNj 

Units 


W (in 2 

st ym) 

(W nr 2 si |irii> DN 

W (m 2 si inn) 

rnsmu 

1 

0 . 452 - 0.514 

0.483 

- 6.2 

191.6 

0 . 77B740 

- 6.98 

1997 

2 

0 519 - 0 601 

0.560 

-64 

196.5 

0.798819 

- 7.20 

1812 

3 

0.631 - 0.692 

0.662 

- 5.0 

152.9 

0.621654 

- 5.62 

1533 

4 

0.772 - 0.898 

0835 

- 5.1 

157.4 

0.639764 

- 5.74 

1039 

5 

1 .547 - 1 748 

1.648 

- 1.0 

31.06 

0.126220 

- 1.13 

230.8 

6 

10.31 - 12.36 

1 1 .335 

3.2 

12.65 

0.037205 

3.16 

N/A 

7 

2 . 065 - 2.346 

2.206 

- 0.35 

10.80 

0.043898 

- 0.39 

84.90 

PAN 

0 515-0 896 

0.706 

- 4.7 

158.3 

0.641732 

- 5.34 

1362 
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Table 5 

ALI spectral range, post-calibration dynamic ranges, and mean exoatmospheric solar irradiance 
(ESUN a ). All EO-1 ALI standard Level 1 products are processed through the EO-1 Product 
Generation System (EPGS). 


EO-1 ALI Sensor (G calmjn = 1 and G c , tmax = 32767) 

Band 

Spectral 

Ramie 

Center 

WaveleiHith 

LMINa 

LMAXji 

^rescale 

^rescale 

ESUNi 

Units 

pm 

W fin 2 st |im| 

|W m 2 sr |im) DN 


W (m 2 4 im) 
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-2.18 
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0.024 

-2.2 

1724 
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0.433 -0 453 

0.443 

-3.36 
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0.045 

-3.4 
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1 
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-4.36 

1405 

0.043 

-4.4 

1996 

2 
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•Katorii 
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0.028 

-1.9 

1807 

3 

0.633 - 0.690 

0.662 

mwm 
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0.018 

-1.3 

1536 

4 

0775 - 0.805 

0.790 

-0.84 
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0.011 

-0.85 

1145 

4P 

0.845-0.890 


-0.641 


0.0091 

-0.65 


5P 

1.200-1.300 


-1.29 

WSEBk 

0.0083 

-1.3 

■Em i 

5 

1.550 - 1.750 

1.650 

-0.597 

31.14 

0.0028 

-0.6 

235.1 

7 

2.080-2.350 

2.215 

-0.209 

29.61 

0.00091 

-0.21 

82.38 
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Table 6 

Earth-Sun distance (d) in astronomical units for Day of the Year (DOY) 
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Table 7 

TM and ETM+ thermal band calibration constants 


Co 11 st tint 

K1 

K2 

Units 

W til * 2 si |jiii) 

Kelvin 

L4 TM 

671.62 

1284.30 

L5 TM 

607.76 

1260.56 

L7 ETM+ 

666.09 

1282.71 
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Appendix 

Table A1 

To maintain consistency, all Landsat scenes are based on the following naming convention 


Format Example: 

LXSPPPRRRYYYYDDDGSIVV 

1 | 
j 

L = Landsat 
X = Sensor 
S = Satellite 

PPP = Worldwide Reference System (WRS) 
Path 

RRR = WRS Row 
I YYYY = Year 

l 

| 

Sensor Examples: 

LM 101 70391 976031 AAA01 (MSS) 

LT 40 1 70361 982320XXX08 (TM) 
LE70160392004262EDC02 (ETM+) 

DDD = Julian Day of Year 
GSI = Ground Station Identifier ' 


W = Version 


j Ground Stations Identifiers - Data received at these sites are held at EROS 

AAA = North American site unknown 

GNC = Gatineau, Canada 

i ASA = Alice Springs, Australia 

LGS = EROS, SD, USA, Landsat 5 data 

1 FUI = Fucino, Italy (Historical) 

acquired by EROS beginning July 1 , 2001 j 

! GLC = Gilmore Creek, AK, US 

MOR = Moscow, Russia 

| HOA = Hobart, Australia 

MLK = Malinda, Kenya 

! KIS = Kiruna, Sweden 

IKR = Irkutsk, Russia 

i MTI = Matera, Italy 

CHM = Chetumal, Mexico 

j EDC = Receiving site unknown 

XXO = Receiving site unknown 

PAC = Prince Albert, Canada 

XXX = Receiving site unknown 


Table A2 

Standard Level 1 product specifications 


Product Type - Level IT (Terrain Corrected) 

Pixel Size - 1 5/30/60 meters 

Output format - GeoTIFF 

Resampling Method - Cubic Convolution (CC) 

Map Projection - Universal Transverse Mercator (UTM) 
Polar Stereographic for Antarctica 
j Image Orientation - Map (North Up) 

Distribution - File Transfer Protocol (FTP) Download only 
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Fig. 1. Comparison of the solar reflective bands RSR profiles of LI-5 MSS sensors. 
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Fig. 2. Comparison of the solar reflective bands RSR profiles of L4 TM, L5 TM, L7 ETM+, and 
EO-1 ALI sensors. 
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